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Abstract 

We seek to derive the probability-expressed in terms of the Hilbert-Schmidt (Euchdean or flat) 
metric-that a generic (nine-dimensional) real two-qubit system is separable, by implementing the 
well-known Peres-Horodecki test on the partial transposes (PT's) of the associated 4x4 density 
matrices (p). But the full implementation of the test-requiring that the determinant of the PT be 
nonnegative for separability to hold-appears to be, at least presently, computationally intractable. 
So, we have previously implemented-using the auxiliary concept of a diagonal-entry-parameterized 
separability function (DESF)-the weaker implied test of nonnegativity of the six 2x2 principal 
minors of the PT. This yielded an exact upper bound on the separability probability of ~ 
0.76854. Here, we piece together (reflection-symmetric) results obtained by requiring that each of 
the four 3x3 principal minors of the PT, in turn, be nonnegative, giving an improved/reduced 
upper bound of § ?s 0.628571. Then, we conclude that a still further improved upper bound 
of « 0.537619 can be found by similarly piecing together the (reflection-symmetric) results 
of enforcing the simultaneous nonnegativity of certain pairs of the four 3x3 principal minors. 
Numerical simulations-as opposed to exact symbolic calculations-indicate, on the other hand, that 
the true probability is certainly less than i. Our analyses lead us to suggest a possible form for 
the true DESF, yielding a separability probability of || « 0.453125, while the absolute separability 
probability of ^^^^^a/a^"^^ ^ 0.0348338 provides the best exact lower bound established so far. In 
deriving our improved upper bounds, we rely repeatedly upon the use of certain integrals over 
cubes that arise. Finally, we apply an independence assumption to a pair of DESF's that comes 
close to reproducing our numerical estimate of the true separability function. 
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Zyczkowski, Horodecki, Sanpera and Lewenstein, in a much-cited article have given 
"philosophical", "practical" and "physical" reasons for studying "separability probabilities". 
We have examined the associated problems which arise, using the volume elements of several 
metrics of interest as measures on the quantum states, in various numerical and theoretical 
studies I2-9I. 

In these regards, we begin our presentation by directing the reader's attention to 
Fig. [TJ The^se depicts various forms of "diagonal-entry-parameterized separability functions" 
as opposed to "eigenvalue-parameterized separability functions (ESF's) 



fD 

8, 



. [TJ These depi 

5SFV) 0, m- 



11 



12[-that we will employ here to obtain estimates and simple exact upper bounds 
on the Hilbert-Schmidt (HS) probability that a generic (nine-dimensional) real two-qubit 
system is separable. 

The subordinate of the three curves in Fig. [l]-derived using an extensive quasi-Monte 
Carlo (Tezuka-Faure [is], Q|) six-dimensional numerical integration procedure-provides an 
estimate of the true, but so-far not exactly-determined DESF. The dominant of the three 



curves-readily obtainable from results reported in 10|, sec. VII]-has the form 
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The intermediate of the three curves, which we first report here, has the same-differing only 
in constants-functional form 

■IJeS (7e^« - 27) t; < 

With each of these three curves we can obtain an associated estimate or upper bound on 
the desired HS separability probability (P^p/reai)- This is accomplished by integrating over 
^ G [—00,00] the product of the corresponding curve with the function (Fig. [2]) (based on 
the jacobian of a coordinate transformation, to be described below) 

_ 64csch^(O(-160sinh(2O - 25sinh(40 + 12^(16 cosh(20 + cosh(40 + 18)) 



that is, 
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Proceeding, thusly, we obtain an upper bound on the HS separability probability of 
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FIG. 1: Three forms of diagonal-entry-parameterized separability functions (DESF's) 
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FIG. 2: Jacobian ([3]), which when multiplied by a separability function and integrated over ^ G 
[— oo, oo] yields the associated Hilbert-Schmidt separability probability 



~ 0.76854 based on the dominant of the three curves, || ~ 0.628571 using the in- 
termediate curve, and an estimate of 0.4528427 for the true probabihty with the subordi- 
nate, numerically- derived curve. (From our work in [12, eq. (25)], we already know that 
the HS probability of a generic real two-qubit system being absolutely separable-that is not 
entanglable by any unitary transformation-is ^^'^^~g'/2^^^ ~ 0.0348338, which then serves as 
a lower bound on the corresponding HS [absolute plus nonabsolute] separability probability 



itself (cf. 
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eq. (29)]).) 

The variable ^ used in the above presentation is the logarithm of the square root of the 
ratio of the product of the 11- and 44-entry of the associated real 4x4 density matrix (p) 
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to the product of the 22- and 33-entries, that is 



^ = logJ = 77 log • (5) 

P22P33 ^ P22P33 



(In our previous studies [Tl, [lO|, we have employed the alternative variables, u = and 



/i = P^P^-^ ^ but now switch to the [more symmetric] form (|5]). Importantly, only the "cross- 

W P'2'2 p33 

product ratio" of diagonal entries is needed in our parameterization to test for separability, 
and not the individual entries themselves.) The jacobian ([3]) used in our calculations is 
obtained by the transformation of one of the diagonal entries, say, P33, to ^ and integrating 
the Hilbert-Schmidt (Lebesgue) volume element (of course, P44 = 1 — pn — P22 — P33) [la, P- 
13646] 

3/3 

dVnS = (PllP22P33P44)~dpiidp22dp33, /3 = 1 (6) 

over pii and P22 and normalizing the result. (To obtain the corresponding HS volume 
elements for the complex 4 x 4 density .natdces. one ,nnst emp.oy-confonning to a p^tern 

familiar from random matrix theory-/? = 2, and (3 = 4 in the quaternionic case (cf. [17]).) 



The use of the celebrated Peres-Horodecki separability test 18|, ll9| is central to our 



analyses. Ideally, we would be able to require that the determinant of the partial transpose 



of p be nonnegative to guarantee separability [20|, |2l|. However, this has so far proved 
to be too computationally demanding a (fourth-degree, high-dimensional) task for us to 
enforce (cf. pij, eq. (7)]). But, in Q, we did succeed in implementing the weaker implied 
test that all the six 2 x 2 principal minors of the partial transpose of p be nonnegative, 
giving us the dominant curve in Fig. [H (Actually, only two of the minors differ nontrivially 
from the analogous set of [nonnegative, of course] minors of p itself.) To derive the sharper 
intermediate curve here, we extended this approach to the four 3x3 principal minors. 
Actually, we found that requiring each of four minors, in turn, to be nonnegative, yielded 
two pairs of identical results. Further, one of these results 

f !^!!:!!(!!!!!z!) .^n 
c (n = } (7) 
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could be obtained from the other set by the transformation ^ — )■ — ^. This curve ([7]) and 
its reflection around ^ = are shown in Fig. 121 The intermediate curve ([2]) in Fig. [H first 
reported here, was constructed by joining the sharper segments of these two curves over the 
two half-axes. A parallel strategy had been pursued with the 2x2 minors. The comparable 
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FIG. 3: The two distinct (red and blue) separability functions obtained from the four 3x3 principal 
minors, the "envelope" (lesser branches) of which defines the intermediate curve in Fig. [TJ The 
y-intercepts of the two curves, are identically ~ 0.867446. 
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FIG. 4: The two distinct (red and blue) separability functions (discontinuous at ^ = 0) obtained 
from the six 2x2 principal minors, the "envelope" (lesser branches) of which defines the dominant 
curve in Fig. [1] 
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results to (171) and Fig. [3] for the 2x2 minors investigation [7] are 

fe-2«(2sinh(0 + cosh(0) C>0 
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and Fig. IH 

For the intermediate curve in Fig. [1] we have the nontrivial y-axis intercept of ~ 
0.867446 (the intercept for the dominant curve being simply 1), while the estimate of the 
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true intercept using the numerically-generated curve is 0.612243, quite close to our previously 



conjectured value of 



1357r^ 
2176 



0.612315 



In obtaining our several results, we used the "Bloore/correlation" parameterization of 
density matrices 22l-l24l| and accompanying ranges of integration-generated by the cylin- 
drical algebraic decomposition procedure 25|, l26|], implementing the requirement that p be 



nonnegative definite-presented in [10|, eqs. (3)-(5)]. The computational tractability of utiliz- 
ing the 3x3 principal minors of the partial transpose in this coordinate frame appeared to 
stem from the fact that each of these four quantities only contains three of the six off-diagonal 
variables (zij) employed in the full parameterization (each set of three variables, additionally 
and conveniently, sharing a common row/column subscript). (The nine-dimensional convex 



set of real two-qubit density matrices is parameterized by six off-diagonal- z 



Pij 



-and 



three diagonal variables-pjj.) Integrating out the three variables not present in the con- 
straint simply leaves us with a constrained (boolean) integration over the cube [—1,1]^, as 
indicated in eq. (3)]-see f|T5|) also. We appropriately permuted the subscripts in the 
indicated coordinate system, so that we could study all four of the minors (thus, finding 
that they fell into two equal sets). Of course, such a simplifying integration strategy is 
not available for the determinant of the partial transpose itself, which contains all the six 
off-diagonal variables (%), rather than simply three. 

Each of the constrained integrations we had initially used, employed as its constraint 
the nonnegativity of a single 2 x 2 or 3 x 3 principal minor of the partial transpose of p. 
(However, we were able above to splice together results, taking the sharper/tighter bounds 



over the half- axes provided by individual outcomes.) We had initial 
either the (Bloore 22|) density-matrix parameterization presented in 



y been unable-using 



10| or the interesting 



partial-correlation parameterization indicated in [23|]-to perform constrained integrations in 
which two or more 2 x 2 or 3 x 3 minors (and a fortiori the determinant) are required to be 
simultaneously nonnegative. (It, then, remained an open question whether or not being able 
to do so would simply lead to the dominant and intermediate curves given already in Fig. [T] 
and by ([T]) and ([2]). However, we were able eventually to answer this question positively for 
the 2x2 minors.) 

We can, however, rather convincingly-but in a somewhat heuristic manner-reduce the 
derived upper bound on the HS separability probability of generic real two-qubit systems 



from II 
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0.628571 to 0.576219 by using a new curve-having a y-intercept of 



- 512 / 



7 



0.752462 as a DESF. (We apply a similar independence ansatz at the very end of the paper 
with quite interesting results [Fig. [9].) This curve is obtained by taking the product of the 
two curves displayed in Fig. [3] (that is, the product of the function with its reflection 
about ^ = 0). A plot of the result shows that it is both subordinate to the intermediate 
curve in Fig.Hl as is obvious it must be, but also clearly dominates the numerically-generated 
curve there, which is an estimate of the true DESF. (Since each of the two curves in Fig. H] is 
simply unity over a half-axis, a parallel strategy in the 2x2 minors analysis can, of course, 
yield no nontrivial upper-bound reduction from ~ 0.76854.) 

The "twofold-ratio" theorem of Szarek, Bengtsson and Zyczkowski j23]-motivated by 
the numerical results reported in j^l-allows us to immediately obtain exact upper bounds, 
as well, on the HS separability probability for generic (eight-dimensional) real minimally- 
degenerate real two-qubit systems (boundary states having a single eigenvalue zero). These 
upper bounds would, then, be one-half those applicable to the nondegenerate case-that is, 
^ 0.38427 and ^ 0.314286. Further, we can, using the results of our numerical 
study, similarly obtain an induced estimate, 0.226421, of the true probability. 

The two sets of derived functions ([1]) and ([2]), based respectively on the 2x2 and 3x3 
minors have the same functional forms, but with differing sets of constants ({1,2,3, 1} vs. 
{9,2048 = 2^^,27,7}). It seems natural, then, to conjecture that the true separability 
unction-which must be based on the determinant of the partial transpose Q, eq. (7)] 
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21| . that is, the single 4x4 minor-will also adhere to the same functional form, but 
with a different set of constants. 

In fact, pursuing this line of thought, as an exercise, we have found that the function 

{315e-3«(-5+18e2«)7r2 
Stb — ^ < 

fits (Fig. [5]) the numerically-generated subordinate curve in Fig. [1] quite well, yielding an 
HS separability probability of f = |l ^ 0.453125, and a y-intercept of ^ 0.6167. 



(Then, by the twofold-ratio theorem 27|, the HS separability probability of the minimally- 



degenerate (boundary) states would be ^ = 2? ~ 0.226563. Also, we have been able to find 
a number of other curves, adhering to this same general structure, fitting the subordinate 
curve in Fig. [1] equally as well, and again yielding || as a separability probability, in addition 
to well-fitting curves yielding somewhat less simple fractions-such as ^ ~ 0.452778, ~ 




FIG. 5: The difference between the numericahy-generated subordinate function in Fig. [5] and a 
hypothetically true separabihty function (l9])-fitting a general pattern observed-giving a separability 



probability of || 0.453125 



0.453086 and |f| ^ 0.45291.) We are obligated, however, to note that in [7|, sec. IX. A] we 



had advanced-based on somewhat different considerations (scaling constants, in particular) 
than her 
equal to 



than here-the hypothesis that this probability is ^ ~ 0.470588, with an associated DESF 



S, 



previous 



(0 



135e-3f (-l+3e2«)7r2 
28-17 
135e«(-3+e2f) 

28TT7 



(10) 



(However, our best numerical estimate at that point was 0.4538838 |lOl . sec. V.A.2] [3, 
sec. IX. A], rather close to our current-study estimate of 0.4528427. By computing standard 
errors of the mean, we can establish a [~ 95%] confidence range-four standard deviations 
wide-for this latter estimate of (0.451634, 0.454051)-that does contain || ^ 0.453125. A 
comparable plot (Fig. |6]) to Fig. [5] shows (|T0|) to provide a considerably poorer fit.) 

One might further speculate-in line with random matrix theory and our previous analyses 
[3]-that the DESF for the generic (15-dimensional) complex two-qubit systems is propor- 
tional to the square of ([9]). If the constant of proportionality were simply taken to equal 
unity, the associated HS separability probability, using the measure ([6]) with /? = 2, would 

306605257r4 



35-52-72-1037r* 



0.252864, rather close to the value ^ ^ 0.242424 conjectured. 



11811160064 230-11 

for a number of reasons, in [7|, sec. IX. B]. Proceeding similarly, using the fourth power of 
iQ, rather than the square, and the measure ([6]) with /? = 4, we obtain the HS quaternionic 

489392789 1 7551757r8 _ 3^° -5^ -7^ -IS-WnSir^ 



probability analogue of 



535315866107766636544 



256-17-19-23 



0.0867454. 
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FIG. 6: The difference between tlie numerically-generated subordinate function in Fig. [5] and the 
previously conjectured true separability function p^ . giving a separability probability of ^ ?a 
0.470588, and a poorer fit than Fig. [5] 

Duplicating the line of analysis of the immediately preceding paragraph, but now using 
the intermediate curve ([2]) given in Fig. [1], instead of the conjectured curve Q-and taking 
the constant of proportionality again to equal 1-we obtain tentative (induced) exact upper 
bounds on the HS separability probability for the complex two-qubit states of I^qqIq^qq ~ 
0.488855 and JotZf^lfiS^foo ~ 0.327414 for the quaternionic two-qubit states. 

In [ij, we studied several two-qubit real, complex and mixed scenarios, in which-in order 
to obtain exact HS separability probabilities-certain of the off-diagonal entries were a priori 
set to zero. In one such (7-dimensional) scenario, we nullified four of the off-diagonal entries, 
allowing only the (1,4)- and (2,3)-entries (the ones interchanged under partial transposition) 
to be complex 0, sec. II. B. 3]. The associated HS separability probability was |. We have 
now been able-parameterizing the off-diagonal entries using polar coordinates-to extend 
this 7-dimensional scenario to a 9-dimensional one, allowing, additionally, any single one of 
the remaining four off-diagonal entries ((1,2), (1,3), (2,4) or (3,4)) to be arbitrary complex. 
The associated DESF is 



^(0 



^e-4«(_i + 4e2e) ^>o 
_ie25 (-4 + e2«) ^ < 



(11) 



with an accompanying HS separability probability of 1| ~ 0.485714. 



Let us now present an additional figure (Fig. [7]) showing-as in Fig. [T]-three DESF's. The 
subordinate curve in this new figure is identically the same as the subordinate numerically 
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FIG. 7: The subordinate curve is the same numerically-derived DESF estimate of the true separa- 
bility probability displayed as the subordinate curve in Fig. [TJ The other two curves are obtained 
by simultaneously enforcing the nonnegativity of certain pairs of the four 3x3 principal minors of 
the partial transpose of p. Both of these superior curves intercept the ^-axis at ^I^^^qq ~ 0.766037. 



estimated curve in Fig. [T] The dominant of the three curves has the form 

7re-fi« (^A/e2«-l(l696e2«-7665e4«-5346ef^« + 188)+3e4«(-7273e2«+1782e'*«+1782) csc-i(e«)) 



71680 



7re-2« (e«-\/l-e2«(~7665e2f+1696e''«+188e6f -5346) +3(-7273e2«+1782e4f +1782) sin-i(e«)) 

71680 ^ < 



(12) 



and the intermediate curve (obtained, as earUer [cf. Fig. |3] by splicing together the ^ < 
and ^ > /esser branches of two curves equal under reflection around ^ = 0), the form 

37r2e-3f ( 18873e2« -4037) 

■ f > u 

(13) 



573440 
37r2e«('l8873-4037e2e^ 



573440 



(This does possess the same basic functional form as already encountered in ([T]), ([2]) and (Q). 
Again, as in Fig. [H of course, the numerically-generated curve yields a separability proba- 
bility estimate of 0.4528427, while the dominant curve yields 0.585542, and the intermediate 
curve, an exact value of = ^^^^ 52 ^ 0.537619. (Both of these last two curves inter- 
cept the ^-axis at ^43^360' ~ 0.766037.) Thus, < || provides a further improved exact 
upper bound on the true Hilbert-Schmidt separability probability. The two superior curves 
in Fig. [7] are obtained by enforcing simultaneously the nonnegativity of certain pairs of the 
four 3x3 principal minors of the partial transpose of p. The dominant curve f[T^ is derived 
by pairing the first and second minors (or, equivalently, three other possible pairs), while 
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the intermediate curve (fT3i) is achieved uniquely by couphng (taking the lesser branches) 
the results pairing the second with the third minor with the (reflection-symmetric) results 
pairing the flrst with the fourth minor. (By the k-th minor we mean the one obtained by 
elimination from the partial transpose of p of its k-th row and column.) 

Given that the direct/naive enforcement of the simultaneous nonnegativity of pairs of 
3x3 principal minors (requiring, a constrained /iwe-dimensional integration) appeared to 
be intractable, we resorted to an alternative strategy to obtain the two superior curves in 
Fig. [71 We exploited the fact already noted above that each of the four 3x3 principal minors 
is parameterized by only three (of the six) off-diagonal Bloore (correlation) variables Zij^s, 
with all of the three sharing a common row/column index (such as the common i index in 
Zij, Zik, Zii etc.), for example, the 4-th minor (with i = 1, j = 2, k = 3,1 = A) takes the form 

minorgxs = 2ehijZikZii - - z^^ - e^^z^ + 1. (14) 

We can, then, arrange-by using a suitably chosen cylindrical algebraic decomposition-that 
any such set of three variables (sharing a common index) comprises the last three to be 
integrated over of the six variables. By performing the flrst (unconstrained) three of the six 
one- dimensional integrations (over, say, Zjk,Zji and zm, in our example), we are simply left 
with (cubical) integrations of the form 

l\ j\ //i^'^^ " '^^^^^ - '^'^^^ ~ zl)dz.,dz.,dza = 1. (15) 

(By reparameterizing the Zi/s in terms of partial correlations [23j, one could re-express 
the full six-dimensional integration as the integral of a simple product measure over a szx- 
dimensional hypercube. But, then, the nonnegativity requirements on the partial transpose 
appear to take on quite cumbersome forms.) Only, at this stage of integration-after having 
integrated out three (extraneous) variables-do we then need to impose (inside the integral 
signs) the three-dimensional nonnegativity requirement of a single minor (minorsxs > 0), 
such as f lT^ to obtain the results reported earlier here. 

To further proceed, in our scheme, we perform the outer two (over zn and Zik, in our 
example) of the three indicated integration steps in f lTB]) -and its analogues-over the cor- 
responding cubes for each of two paired minors independently of one another, and then 
combine (multiply) the two results together, which are then integrated (in a joint manner) 
over the remaining shared/ last variable {zij in our illustration here) to derive the new curves 
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FIG. 8: The two distinct (red and blue) separability functions obtained from enforcing the joint 
nonnegativity of the first and fourth 3x3 principal minors of the PT-giving the dominant (blue) 
curve on the left-as well as the joint nonnegativity of the second and third 3x3 principal minors 
of the PT. The lesser branches of the two curves define the intermediate curve (|13p in Fig. [71 The 



in Fig. [71 Our approach here, thus, consists in replacing a direct [but intractable] five- 
dimensional constrained integration-five being the number of variables parameterizing any 
two of the four 3x3 principal minors, we want to be simultaneously nonnegative-by a pair of 
independent constrained two-dimensional integrations (each member of the pair concerned 
with the nonnegativity of only a single minor) conducted over three-dimensional cubes. The 
two distinct one-dimensional results obtained are, then, joined by multiplication together 
into a single one- dimensional integration (over Zij, the shared variable, in our illustration). 
Since there is a factor of (|)^ = || in the three- fold integrals f[T^ . we importantly assign-by 
symmetry-a weight of | to each single-fold integration step taken. 

The intermediate curve in Fig. [TJ given by (fT3l) . is constructed by taking the lesser 
branches of the two curves in Fig. [H] (The four possible pairings of minors other than the 
first with the fourth, and the second with the third all yield the same dominant curve shown 
in Fig. [7l) The two greater branches in Fig. [H together yielding an upper bound on the 
separability probability of — ^l^i'f ~ 0.854936, take the form 



573440 S ^ " 

Let us-similarly as we have done before-take as a new separability function the product 



greater branches are described by (I16p . The value at the intersection is 



111277r^ 
143360 



0.766037. 




(16) 



13 



Sep. funct. 





0^ 
















/ 0.4 






1 , , // , 


// 0.3 






-2 -1 


1 ^ 





FIG. 9: The (blue) separability function derived by taking the product of the two curves displayed 
in Fig. [8] along with the numerical estimate of the true separability function. The latter (red) 
curve crosses the y-intercept at 0.612243 and the (blue) product DESF at the lesser value of 
0.586813. 



1238101297r"t 
2^-4.52.72 



of the two curves displayed in Fig. [HI Use of this product DESF in formula (jl]) yields 



7f 



(1803179l7r2 - 177044420) 



Kiireai = — . 52 . 72 ' ^ 0.453503, (17) 

very close to our earlier numerical estimates of 0.4538838 10, sec. V.A.2] [t!, sec. IX. A] and 
lying within the confidence range (0.451634, 0.454051), we established above. (Possibly, this 
product DESF is, in fact, the function that would arise if one could simultaneously enforce 
the nonnegativity of a// four 3x3 principal minors [but see final paragraph]. However, it lacks 
the simple functional form repeatedly previously observed above.) We display this derived 
product separability function in Fig. IHlalong with the closely-fitting numerical estimate of the 
true function, already appearing as the subordinate function in both Figs. [H and [71 (If we, 
in a similar vein, take as a product DESF, the square of f[T2|) . that is, the dominant curve 
in Fig. [3-since this arises identically from four of the six possible pairings of minors-the 
associated separability probability falls, rather unrealistically to 0.367762. So, independence 
of minor pairings does not appear to be a tenable hypothesis in this case.) 

It does, however, appear that we can reduce the y-intercept in Fig. [7] from ^H^^qq ~ 
0.766037 to if|§ ^ 0.6872 by enforcing the simultaneous nonnegativity of the second, third 
and fourth minors, using repeated integration over cubes. 
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